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Abstract. We study the Cauehy problem for one-dimensional dispersive system of Boussi- 



1. Introduction 



S . nesq type which models weakly nonlinear long wave surface waves. We establish the local 

well-posedness and ill-posedness of solutions to the system. We also provide criteria for 
the formation of singularities. 

(N 

< 

^ ■ The theory of water waves embodies the equations of fluid mechanics, the concepts of 

wave propagation, and the critically important role of boundary dynamics. Due to the 
complexity of the governing equations for water waves, physicists and mathematicians are 
led to derive simpler sets of equations likely to describe the dynamics of the water waves 

^ ■ in some specific physical regimes. 

\^ • One regime that arises in practical situations is that of waves in a channel of approx- 

^ . imately constant depth h that are uniform across the channel, and which are of small 

amplitude and long wavelength. Similar waves appear as long-crested disturbances on 



larger bodies of water. Let a be a typical wave amplitude and A a typical wavelength, the 



o 

CN ' conditions above amount to 

X- 

^ , The equations within the above scaling regime couple the free surface elevation r] to 

- - -' the horizontal component of the velocity v. When the nonlinear and dispersive effects 

are balanced, that is, e = 0{6), and for one dimensional surfaces and fiat bottoms, at 

the asymptotic expansion to the first order in e, the system reduces to the well-known 

Boussinesq equation [HI [2T1 1231 EI] 

(note that rj = v + 0{e)) and any perturbation of the surface splits up into two components 
moving in opposite directions. 
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2 R. M. CHEN AND Y. LIU 

Taking advantage of the freedom associated with the choice of the velocity variable and 
making full use of the lower-order relations (the wave equation written as a coupled system) 
in the dispersive terms, Bona- Chen- Saut [9] put forward the following three-parameter 
family of Boussinesq systems (referred to as the ahcd system) for one dimensional surface 
(a two-dimensional analogue is derived by Bona-Colin-Lannes [H]), 

^Vt + Vx + {'nv)x + avxxx - br]txx = 
\ ) 



Vt + Vx + vvx + crjxxx - dvtxx = 0, 

which is formally equivalent models of solutions of the Euler equations. In the above system, 
ri{t, x) is proportional to the deviation of the free surface from its rest position, and v{t, x) 
is proportional to the horizontal velocity taken at the scaled height 6 with 0<6<1{6 = 1 
at the free surface and 6* = at the bottom), and 

with V and /i arbitrary real numbers, and r > the surface tension. Throughout this paper 
we will just be dealing with the pure gravity waves, that is, r = 0. 

Preliminary work on well-posedness of (II. ip is detailed in [8] for various collections of 
a, 6, c, and d. The authors proved that all the linearly well-posed systems are locally 
nonlinearly well-posed. Notice that an interesting feature of (11.11) is that its dispersive 
properties can be quite different according to the choice of a, 6, c, and d. In fact the 
dispersive relation is 

,.2..^_ |.|2 (i-«ien(i-cien 

and hence the (a, c) pair enhances dispersion while the (6, d) pair weakens it. In this paper 
we will concentrate on the case of weak dispersion, that is, a = c = and b = d = |. After 
rescaling b and d to 1 the system becomes 

(1.2) <^ ^ X eR, t>0. 

[ Vt - Vtxx + Vx +VVx=0, 

One of our goals is to obtain a sharp well-posedness result of the above system. It 
has been shown in [^ through a contraction principle that system (II. 2p is well-posed in 
Sobolev spaces i7*(M) x H'^lM.) for s > 0. Using the idea of Bejenaru-Tao [3] we manage 
to prove, via explicitly constructed initial data, that the flow map fails to be continuous 
below L^(R) X L^(M). The idea stems from looking at the second order iterates in a Picard 
iteration about the chosen initial data, which leads to a high-to-low frequency interaction: 
a solution starts off initially with small energy and Fourier transform supported primarily 
at high frequencies, but quickly generates a large energy at low frequencies. This norm 
explosion indicates that the solution operator, if it exists at all, has a severe singularity 
near the zero solution. 
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Our ill-posedness threshold agrees with the one for the BBM equation [10]. However, 
system (II. 2p involves both the surface rj and velocity v components and consequently gen- 
erates many more interactions in the second iteration map than in the scalar BBM case. 
Our initial data is chosen to minimize the surf ace- velocity interactions but to maximize the 
frequency interactions of the linear propagators in the nonlinear terms. Similar techniques 
were used in proving a C^ ill-posedness of the full gravity-capillary water wave problem 

m- 

Note that the weak dispersive system (11. ip corresponding to a = c = does possess 
global-in-time solutions. For the particular value ^^ = |, the system features exact solitary- 
wave solutions [13]. However for general 6 G [0, 1] the existence of solitary waves is still 
open. 

From modeling point of view, it is interesting to consider the following scaled asymptotic 
model of (OD 

Vt + Vx + e [(r]v)x - ijtxx] = 0, 



'1-3) , 

[vt + r]x + e [vvx - vtxxl = 0. 

A natural but important question about the asymptotic water wave models is the jus- 
tification of their hydrodynamic relevance, which involves two issues, whether or not (i) 
solutions to the full water wave equations which satisfy the modulational Ansatz exist for 
a sufficiently long time, and (ii) solutions to the asymptotic models approach the solutions 
to the full water waves equations within certain time. 

It was proved in [2] that the existence time of the full water wave equations in the long 
wave regime {e = 0{6)) is order 0{l/e), and hence the large-time existence of solutions 
is provided. In regard to the second question, it was shown in [9], [13] that all Boussinesq 
models give an approximation of order 0{eH) to the full water waves equations, provided 
the solutions exist. Therefore, the justification of system fll.3p requires a large time well- 
posedness of the Cauchy problem. 

For two dimensional surfaces, it was shown by Bona-Colin-Lannes [H] using a hyperbolic 
approach that the existence time of solutions to system (11. ip is order 0{l/e) in H2 + (M.'^) 
when a = c and the system is symmetric. In the strongly dispersive case when a = c = 1/6, 
with b = d = 0, and r = 0, it was proved by Linares-Pilod-Saut [25] that the solutions to 
system (II. ip exist in a time scale of order 0{l/^/e) in H2~^(M.'^). On the other hand, it 
was also proved in [22] that when b,d > system (II. ip is well-posed in the Sobolev space 
H^{M.'^) on the time scale of order 0(£~2+). The only success in reaching the optimal time 
scale for general a, b, c and d, to our knowledge, is due to Ming-Saut-Zhang [28] for a, c < 
and b,d > 0, where a Nash-Moser approach is applied when the initial data are smooth 
(and with a loss of derivatives). 

For one dimensional surfaces, for properly chosen initial data t]q and vq with vq = rjo + 
0{e)riQ, it was proved in [T] that the solutions to system (II. 2p persists on the time scale of 
order 0{l/e) with an extra high regularity assumption. One of our goals here is to obtain 
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the well-posedness of system (11. 3p {a = c = 0, b = d > 0) on time scales of order 0(1 /y/e) 
for general initial data in a much rougher space, achieving the rigorous justification of 
the model on the corresponding time scales. In general, we, however, are not able in our 
functional setting to reach the optimal time scales of order 0{l/e). 

Another aspect of this paper concerns the formation of singularities of system (11. 2p . 
Applying an energy-type argument we show that solutions blow up in finite time when v^ 
is not bounded below. Together with the use of the Hamiltonian functional, we further 
prove that blow-up in finite time happens if and only if t] becomes unbounded below, which 
improves the result of global existence obtained in [1] and also seems to agree with some 
numerical predictions in |15] . 

Since dissipation effects are ignored in the derivation of (II. ip and the overlying Euler 
system is Hamiltonian, it is expected that some special cases of (II. ip will also possess a 
Hamiltonian form. In particular, for our case system (II. 2p . the following functional 

(1.4) niv, v) = l- f [r]^ + {l + r]y] dx 

serves as a Hamiltonian. In addition, system (II. 2p has the following conserved quantities 

T] dx, / V dx, / {rjv + rj^Vx) dx 



along with 7/(77, t>) ([S])- 

The plan of the paper is as follows. In Section 2, we introduce the notations and function 
spaces to work with. In Section 3, we establish the analytical well-posedness of system (II. 2p 
in Sobolev spaces H^{M.) x if'*(M) for s > 0, and we further show that the time scale of 
existence is of order 6~^'^. In Section 4, we prove that system (II. 2p is C^ ill-posed below 
L^(R) X L^(M). In Section 5, we discuss the formation of singularity and provide a series 
of blow-up criteria. Finally in Appendix, we give the decay estimates of the linearized 
equation. 

2. Resolution space 

In this section we introduce a few notation and we define our functional framework. 

Let u{t,x) = {r]{t,x),v(t,x)). Denote \u(t)\p^q = |'7(^)|p+ K'('^)l'? the norm in L^ x L'^, 
for 1 < p, q < +00, where | ■ \p refers to the L^ = L^(M, dx) norm. Let ||u(t)||s = 
||?7(t)||s + ||w(t)||s for s > 0, denote the norm of Sobolev space X^ = H'^(R) x H^(R) where 
II M lis is the norm of Sobolev space H^{M.). 

For u = u(t, x) G iS'(M^), let u be its Fourier transform in space. We define the Japanese 
bracket (x) = (1 + Ix^)"*^'^. The notation A r^ B means that there exists a constant c > 1 
such that -\A\ < \B\ < c\A\. For any positive A and B, the notation A < B (resp. A> B) 
means that there exists a positive constant c such that A < cB (resp. A > cB). 
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Notice that one may rewrite ( II. 2p as the following abstract form 

(2.1) ut = Au + N{u) 
where 

(2.2) ^ = -(1 - d'J~'d. (J J) , N{u) = -(1 - dlr'd^ (p) . 
Then the linear semi-group S{-) associated with (II. 2p is given by 

(2.3) (isi) (0 . ( -(i<|p, -"«|>::;)) h). 

In this way, we will mainly work on the following integral formulation of (II. 2p with initial 
data Uo 

u{t) = S{t)uo + j S{t- t')N (n(t')) dt' 

where for u(t) = {ui,U2)(t) and v(t) = (fi,f2)(^) 

(2.5) N,iu,v) = -\f^S{t-m-dird. {^"XvT") (^') ^'' 

is a bilinear operator. 

We can apply a Picard iteration (fixed point argument) to solve the above integral 
equation locally in time in the space X^ = C([0,T],X*) equipped with the usual norm 

ll^lUf, = sup||u(t, ■)||,. 

[0,T] 

Clearly, the free evolution operator S(t) is unitary on X'^ for any t > 0, that is 

\\S{t)uo\\s = \\uo\\s, 
and hence for any T > 0, 

(2.6) \\S{t)uo\\x^ = \\uo\\s- 
Similarly, 

(2.7) \\N2iu,u)\\x^ < T\\Sit)N{u{t)) \\x^ = T\\N {u{t)) |U- 

Therefore the Picard iteration can be used if we can establish a bilinear estimate of 
\\N{u{t))\\x^. 
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3. Local well-posedness 

The goal of this section is to estabhsh the local analytical well-posedness of system (ll.2p . 
As is explained above, we would first need a bilinear estimate for the nonlinear part in the 
integral formulation (12 ■4p . This is achieved using the following lemma , which can be found 
in m. 



Lemma 3.1 ([lO]). Let u,v E H'{R), s > 0. Then 

(3.1) 11(1 - dlr'd,{uv)\\Hs < Cs\\u\\h44h^, 
where Cg > is some constant depending only on s. 

Applying the above lemma to fl2.7p we obtain 

(3.2) \\N2iu,u)\\x^<CsT\\u\\l^. 

The following lemma will be used to estimate the maximal time of existence for solutions. 
Lemma 3.2. Let s > s' > 1/2, s' > s - 1 and u,v e H^{R). Then 

(3.3) Wil - dlr'd,{uv)\\Hs < Cs\\u\\hs'\\v\\hs, 
where Cs > is some constant depending only on s. 

Proof. First it is easy to see that 

\\{1 - dl)-'d.,{uv)\\Hs < \\uv\\hs-u 

When s > 1, from [21] and the assumption that s' > 1/2, s' > s — 1 we know that the 
elements of LL^ (M) are multipliers in H^~^{M.), that is 

and hence we have (13.31) . 
When s < 1, 

||Mt;||j^s-i < ||Mf||^2 < ll^ll^s' 11^11^2 < ||M||/^s'||f 11//", 

thus we also arrive at (13. 3p . D 

3.1. Bejenaru-Tao's abstract well-posedness theory. With the help of (12. 6p and (13. 2p . 
the main analytic well-posedness theorem is given as below, where we follow the idea es- 
tablished by Bejenaru and Tao [3]. For a complete statement and proofs of the abstract 
result please refer to [3]. Here we just provide an argument for the estimate of the maximal 
existence time. 

Theorem 3.1. (Local well-posedness in X^ = L^(M) x L'^{M.)) Fix s > 0, for any iio(x) = 
{riQ{x) , Vq{x)) G X'^ , there exist a T = T{uq) > and a unique solution u{t, x) = {ri(t,x),v(t,x)) E 
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X^ to equation (12. 4p . and hence to (11.21) with initial value Uq. The maximal existence time 
Ts for the solution has the property that 

4:Cs\\Uo\\s 

where Cs > 0, and 

(3.5) Ts < oo ^ limsup \\u\\s = oo. 

Moreover, Tg is independent of s if s > 1/2. 

More specifically, we can construct the iterations An : X"^ — )■ X^ for n = 1,2, . . . by the 
recursive formulae 

Ai{uo) = S(t)uo 

(^■^) An{Uo)= Yl N2{An,{Uo).An,{Uo)) for 71 > I, 

ni,n2>l,ni+n2=n 

with the property that 

\\Anif)\\x^ < K"" ll/llx» ' for some K > 0, 
and we have the absolutely convergent (in X^) power series expansion 



oo 



(3.7) u = J2Muo) 



n=l 



This implies the analyticity of the solution map Uq^^ u from bounded sets of X^ to X^. 

Proof. We only provide an argument for the maximal time of existence here. For any initial 
data Uq G X'^, define a map 

F : D = {u e X^ : ||n||^s < R = 2||no|U4 -^ X^ 

as 

F{u) = S{t)uo + N2{u,u). 
Then from (12.61) and (13.21) . for any u,w E D, 

\\F{u)\\xs^ < \\uo\\s + CsT\\u\\j,^ <^ + CsTR\ 

\\F{u) - F(tU)|Uj < CsT\\u + tv\\xs^\\u - tv\\x^^ < 2i?CJn - w\\x'^. 
Hence by choosing 

T ' ' 



2C;R 4C.||So||. 

the map F becomes a contraction from D into itself. Therefore for each s > 0, the maximal 
existence time Tg satisfies 

T> ' 
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and Tg is finite implies that the X'^-norm of solution blows up to infinity as t — )■ T^. Hence 
we have proved (13.41) and (13. 5p . 

It is obvious that Tg is a non-increasing function of s. Let s > 1/2. For any 1/2 < s' < s, 
suppose that Tg < Ts'. Therefore by uniqueness the maximal X^ and X^ solutions coincide 
on [0,Ts). Also from the definition of T^ and Tg', ||'u(t)||s/ remains bounded on [0,Ts] while 
||'u(t)||s blows up as t — !■ Ts. From (12.41) and Lemma [3^2] we have 

\\u{t)\\s< \\uo\\s+ [ \\N {u{t')) \\s dt' 



(3. 



< 11% I 



diro^ 



'^ \ ly2 ] (f) 



dt' 



< ||^o||. + C, / \\uit')\\s'\\uit')\\sdt'. 



Then by Gronwall's inequality we conclude that ||-u(t)||5 is bounded on [0,Ts), which is a 
contradiction. Therefore, Tg = Tg' and hence it does not depend on s when s > 1/2. D 



3.2. Hydrodynamic relevance. As is discussed in the Introduction, it is of interest from 
modeling point of view to consider system (II. 3p . 

The result of Theorem 13.11 does not indicate the dependence of the temporal interval 
[0,T] on the parameter e. From |2], the physically relevant temporal regime for (II. 3p is 
from 0{l/e) up to 0{l/e'^). Hence we present here the mathematical analysis on that 
perspective. 

In an abstract form, the above system can be rewritten as 



where 

A' = -(1 - ed^r^d:, 

The associated linear semi-group 5'^(-) is 



ut = A'^u + N^u) 
° ^ ' NHu) 



<i - edir'd. ( p 









m 



which is still unitary on X'^. 

The Duhamel formulation of (11.31) is 



u{t) = S%t)uo + / S%t - t')N' {u{t')) dt'. 
Jo 



The bilinear estimate is 
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Lemma 3.3. Let u,v ^ H^ , s > 0. Then 

(3.9) ||£:(1 - edl)~'^d,^{uv)\\Hs < CsV£\\u\\hs\\v\\hs, 

where Cg > depends only on s. 



Proof. Using duality and a polarization argument, one may write (13. 9p in the equivalent 
form 



(3.10) 



e^iO' 



M^iM^-^iMOd^d^i 



< Csy/£\\u\\L2\\v\\L2\\w\\L2. 



;i+^a(ei)^(e-ei)^ 

For s > 0, {0^ / {^i)"^ {C, — ^lY < Cs and hence such term in fl3.10p can be ignored. Let 



MO 



l + eC 



MO- 



Then the left hand side of fl3.1Up is {u * v,wi). Moreover, since u{x) and ui{x) = u{—x) 
have the same L^-norm, we have 

{u * V, wi) = (^ui * wi, v') < \\ui * u;i||^2||'u||^2 
< IklUHI^ilUilK'IU^ 
si 



< 



i + ee 

and hence proves the lemma. 



^||L2|K'llL2|klU2 < v^||n||L2||t;||i2||w;| 



L2, 



L2 



D 



From the above lemma we have for s > that 

ft 



S%t-t')N'{u{t')) dt' 



< C.v^T||n||^.. 



X;; 



Consequently we obtain the existence of solutions on the time scale of order 0{1/ ^/£). 

Theorem 3.2. Fix s > 0, for any uq{x) = {rio{x) , vo{x)) G X**, there exist aT = T{uq) > 
and a unique solution u{t,x) = {r]{t,x),v{t,x)) G X^ to f lL3p with initial data uq. The 
maximal existence time T = Tg for the solution has the property that 



(3.11) 



T, > 



1 



4:^/eCs\\uo\ 



where Cs > depends only on s. When s > 1/2 the maximal time T, is independent of s. 
Moreover, the solution map uq \-^ u is analytic from bounded sets of X^ to X^. 
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4. Ill-posedness 

The ill-posedness result we discuss here can be viewed as an application of a general 
result proved in |3]. Roughly speaking, this general ill-posedness result requires the two 
following ingredients: 

(1) The equation is analytically well-posed until a certain critical index Sc and the 
corresponding solution-map is also analytic in the space of solutions. 

(2) Below this critical index, one iteration of the Picard scheme is not continuous. The 
discontinuity may be driven by high-to-low frequency interactions that blow up in 
frequencies of order smaller or equal to one. 

The first ingredient is given by Theorem 13.11 whereas the second one can be derived using 
the abstract ill-posedness theory in [3] (c.f. Proposion 1 in [3]), which says that if the map 
uq I— !■ J2'^=i ^n(wo) is continuous in a coarse topology, then each component uq H- An{uo) 
of the series is also continuous in this coarse topology. Therefore, we may reduce the 
discontinuity of the flow map to disproving a bilinear estimate, namely, we will show that 
the second iteration A2{uo) is not continuous. 

Our ill-posedness result is stated below. 

Theorem 4.1. (Ill-posedness below X^) For any s < 0, let R> be arbitrary. Denote Br 
the ball 

Br = {no e X' : ||no||. < R}- 
Let T and uq ^ u be as in Theorem \3.1[ Then the solution map uq ^ u is discontinuous 
from Br (with the X^ topology) to Xj, (with the X^ topology) for any s' G M. 

Remark 4.1. Note that for the one-component analogue of (11.21) . c.f. the BBM equation 
[1] , it was shown by Bona-Tzvetkov [TU] that the associated Cauchy problem is C^ ill-posed 
below L^(]R) by proving that the second Picard iteration fails to be continuous. In fact 
their problem fits in the framework of Bejenaru-Tao, and hence one can infer that the 
BBM equation is C^ ill-posed below L^. Our result here suggests the same threshold for 
the system (11.21) . but the frequency interactions here are more delicate. 

4.1. Second iteration A2{uo)- By definition (13. 6p . we know 

A2{uo) = N2{S{t)uo,S{t)uo) = [ S{t-t')N{S{t')uo) dt'. 

Jo 

Our goal is to construct initial data Uq that violates 
(4.1) \\A2{uo)\\x^<\\^ofs ioTsKO. 

For simplification of notation, we denote 

^r]\ _ (LiTj + L2V' 



^ ' ^ V I \L2Ti + LiV 
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where 

(4.2) 



Li(e,t) = cos(A(Ot) 
L2(e,t) = zsin(A(0t), 



with A(.^) = ^ {^) . For an initial datum uq 



A2{uo)= / S{t-t')il 



:;»,let :;; W^SW -l.Then 



Q] 



«-'a. (1:5) (0 *' . / (^: 1 *'^ 



Taking the Fourier Transform of the above second order iterations yields the following set 
of terms 



Qi 



^^ 



1+e 



L^{U-t')rn^i{U') + L2{U-t')^vl{U') 



Li(e,t-ou^i(ei, 0^0(6) + ^2(6, 0^0(^1 



X L2(e-ei,t')^o(e-6)+^i(e-ei,0^o(e-^i)Mei 



^^ 



L2(e, t - t') L2(ei, t')Uii) + ^1(6, t')vo{ii 



2(1+ a 

X ( L2(e - Ci, t')Ui - 6) + L,{i - ^1, t')^o(e - ^1) ) ^6, 



Q2 



^e 



i+e 



22(e,t-t')T^i(e,t') + Li(e,t-to^^?(e,o 



L2(e,t-t') Li(ei, 0^0(6) + ^2(6, t')t^o(ei 



X L2(e-ei,o^o(e-6) + ^i(e-ei,o^o(e-ei)Hei+ 



^e 



Li(e,t-t')U^2(ei, 0^0(6) + ^1(6, 0^0(^1 



2(1+ a 

X ( L2(e - ei, t')Ui - ii) + ^i(e - ei, t')t^o(e - eO ) di^. 



Of the above 16 terms in the two Duhamel operators Qi and Q21 we look for a dominant 
lower bound term that violates an assumed upper bound in 04.11) . This lower bound will 
arise from the following construction. 
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4.2. Choice of data. For a large A^ we choose initial data Uq = (t/c'^o) so that 

^0 = 0, 

where x is the characteristic function. In this way we eliminate all interactions in A2{uo) 
involving rjQ. Then our expressions for Qi and Q2 reduce to the following: 



(4.3) 



Qi 



^^ 



1 + e 



Li(e, t - t')i^2(6, t')L^{^ - iut') + -L2(e, t - f)L^{ii, t')L,{i - 6, t') 
X t^o(6)^o(^ - 6) d^i. 



Qi 



^i 



l + e 



X ^o(6)^o(^ - 6) c?^i- 



Proposition 4.1. //A^ zs sufficiently large, then on the support of Vo{^i)vo{^ — ^i) and for 
<t' <t < 1 we have the following bound: 

L2(e,t- 0^2(6, 02i(e- 6,0 + ^2i(e,t-t')Li(6,t')2i(e- 6,0 >^- 

Proof. Since A^ is large we know that on the support of t'^o(6)^o('C ~ 6) 

^<|A(6)|,|A(e-6)l<^. 



So when < t' < t < 1, 



l>Li(ei,t'),^i(e-6,0> 



As for L2, we know that 

£2(6,0 =|sin(A(6)0l< 



sm 



t' 1 
< — < — . 

- A^ - A^ 



Hence similarly 



i^2(e-6,0 



1 
< — . 

- A^ 



Because on the support oi vo{^i)vo{^ — ^i), N — ^ < |6|, I'C ~ 61 < A^ + |, we deduce that 
16 < 1 or 2A^-1 < 16 < 2A^ + 1. 

If 16 < 1, then 

\m\<\m\<i 
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So 



Li(e, t-t') = cos (A(0(t - t')) > cos ( ^) > ^, 



L2{^,t-t'] 



|sin(A(0(t-tO)|<sin(i]<i. 



Therefore we have 



L2(e,t-t')22(ei,o2i(e-ei,o 



11 _ 1 

- 2 ' iV '^ ~ 2iV' 



2 ^' ' ^' ' ^' '"''-2222 16 

Hence by choosing A^ sufficiently large, say, A^ > 16 we obtain the desired bound. 
If 2A^ - 1 < 1^1 < 2iV + 1, then |A(OI < ^ and hence 



L,{U-t')>^, \L2iU-t'] 



< 



1 

N' 



Therefore 



L2(e, t -t')^2(6, 0^1 (e-ei,o + 2^1 (e,t- 0^1(^1, o^i(e-ei,o 





1 
> 

- 16 


1 1 

2iV2 - 32 


for iV > 4. 







D 

Remark 4.2. Here we only estimate the symbols in Q2 because we are seeking a lower 
bound of 11742(1*0) 1 1 x"'- From the definition of the symbols Li and L2 we see that a lower 
bound can be obtained only when there is a multiplier consisting of all Li-symbols. This 
is why we assign rjQ = 0, since iivo = then the corresponding Qi and Q2 are 



Qi 



^^ 



i + e 



X ^o(^i)^o(^ - 6) d^u 



Q2 



^e 



i + e 



L2(e, t - t')^i(ei,t')^2(e- 6, o + 2^i(e,t- 0^2(^1,0^:2(^-6,0 



X ^o(6)^o(^ - 6) d^i, 
and all the multipliers contain the symbol L2. 
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4.3. Ill-posedness. Choosing initial data as in (14. 3p . we have 



(4.4) 



Fol 



Pol 



H^ 



Ker^oL^-i 



Proof of Theorem 4-i Fix s < and s'. Suppose for contradiction that the corresponding 



solution map is continuous. Then from Proposition 1 in [5], we have (14. ip . 

From (13. 4p and (14. 4p we may rescale T to be 1. Consider a low- frequency region I 
{I'd ^ 1}) cind therefore (^)* ~ 1. It then follows from Proposition 14. 1 1 that 



^2(^0) \x^' = s^P 
^ o<t<r 


Jo KQ^J 


dt' 


qxq 


= sup 

o<t<r 


iO' f Qidt' 
Jo 


+ sup 
q o<t<T 


> sup 

0<t<T 


iO' fQ2dt' 
Jo 


^l 


1 
> 
- 32 


J- + ^ Jo JR 


^N-^ \^\lhi)^N-% 


which violates f4.ip when s < 1 


]. 









Jo 



r^ 
^i 



Lj{l) 



n 



5. Blow-up criteria 

In the previous two sections we establish the local well-posedness and ill-posedness of 
solutions to system (II. 2p . Attention is now given to the formation of singularities. 

As is pointed out in |8], (II. 2p possesses a Hamiltonian structure with the Hamiltonian 
functional 



■^ JR 



and hence for any local solution u G X* with s > 1/2, ^{(n) is conserved. 
Now we give our first blow-up criterion. 

Theorem 5.1. Let s > 3/2. A solution u = {rj^v) of (II. 2p with initial data Uq G X* blows 
up in finite time T < 00 if and only if 



(5.1) 



lim inf 



inf Vxit, x) 



-00. 
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Proof. We first prove the result for some s > 2. Multiplying (ll.2p by m and using integration 
by parts we obtain 

III* = I /.(■"' + ''"' + "' + ■"3''^ 

(5.2) ^-/M + M.„ + ,.. + «V]d. 

= - Vxtf dx. 



Thus applying Gronwall's inequality we know that the X^-norni of the solution does not 
blow-up in finite time if and only if Vx{t, ■) is bounded from below on [0, T). 

Next we use fl3.8p in Theorem 13.11 to conclude that the X^-norm of solution does not 
blow up for all 1 < r < s if and only if Vx{t, ■) is bounded from below on [0, T). Therefore 
we proved the statement of the theorem for all s > 2. Lastly, a density argument and the 
continuous dependence on initial data ensures the validity for all s > 3/2. D 

Remark 5.1. The result in Theorem 15.11 indicates the blow-up of the slope of velocity. But 
it is quite different from the phenomenon of "wave-breaking" (solution remains bounded 
but the slope becomes infinity in finite time) which happens in a class of shallow water 
models [121 [13 [HI [IHl 1201 ISTJ [32]. For our system (11. 2p . the conservation laws are not 
strong enough to provide a priori bounds for the solutions. Though in Theorem 15.11 the 
regularity index is assumed to be above 3/2, from Theorem 13.11 we see that when s > 1/2 
(hence the functions are pointwise defined), the maximal time of existence is independent 
of the regularity index s. Hence this suggests that one may be able to look for a blow-up 
criterion in a lower regularity setting. In fact we can further refine our blow-up criterion 
as below, which asserts that it is indeed the pointwise blow-up of the solution itself rather 
than its slope. 

Theorem 5.2. Let s > 1/2. A solution u = [rj, v) of (ll.2p with initial data uq G X'^ blows 
up in finite time T < oo if and only if 



(5.3) liminf 



i-s>T 



inf riit.x) 

xeM. 



-oo. 



Proof. "^": Similarly, we first prove the result for s > 1. If not, suppose that ?7(t, x) > —M 
on < t < T for some M > 0. It follows from the conservation law of the Hamiltonian "H 
that 

I rf dx = Uiuo) - / (1 + r])ii^ dx < Hiuo) + M v^ dx. 

J«. Jr Jr 
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On the other hand, multiplying the second equation in (11.21) by 2v, integrating over 
and using the above inequality we have 

v"^ + vl dx = 2 j v^r] dx < \vx\l + \ri\l 



<'H{uq) + M\v\1 + \v^\1 
<n{uo) + {M + l)\\v\\l,. 
Then using Gronwall's inequality we deduce 

\\v\\m < C{uo,T). 

Now multiplying the first equation in fll.2p by 2?7, integrating over M, and using the 
estimate for \\v\\h^ we have 

V'^ + Vx dx = —2 I rjVx dx + 2 I r]r]xV dx 



dt _^ 

< '^Ivxhlvh + l^^loo (|^|2 + \Vx\l) 

^ \\v\\h^ + Mm- 
Hence 

\\r]\\m < C{uo,T). 

In this way we see that the X^-norm of the solution does not blow up. Then from (13. 8p 
in Theorem 13.11 we conclude that the X'^-norm of solution does not blow up, which is a 
contradiction. 

Now for 1/2 < s < 1, a density argument and the continuous dependence on initial data 
ensures the same result. 

"<^": This is can be proved by a use of the Sobolev embedding H^(M.) "-^ L°°(M) for 
s>l/2. D 

Remark 5.2. (1) From the above two theorems it is clearly seen that when s > 3/2, the 
blow-up of the slope of the velocity is equivalent to the blow-up of the surface profile. 
(2) From Theorem 15.21 we also obtain a criterion for global existence of solutions: for 
s > 1/2, if ri(t,x) is bounded from below on t G [0,Tmax) , then Tmax = oo. This is 
an improvement of the result in [1], where the global existence of solutions requires that 
r]{t, x) > —1 + a for some a > 0. 

Finally we give a blow-up characterization of the velocity profile. 

Corollary 5.1. Let s > 1/2. A solution u = {ri,v) of (11. 2p with initial data uq G X* blows 
up in finite time T < oo if and only if 

(5.4) limsup |f (t, ■)|oo = oo. 

t-s>T 
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Proof. "^": If not, suppose \v\oo < M on < t < T for some M > 0. From (I5.2p 



d 
~dt 



\u' 



li = — / Vx'rf dx = 2 vrjrjx dx < 2M||?7||^i. 

Jr Jr 



So by Gronwall's inequality, \\u\\i is bounded. Hence |?7|oo remains bounded, which is a 
contradiction. 

"<^" : Again we can prove this by the Sobolev embedding H^(R) M- L°°(M) for s > 1/2. D 



Appendix A. Linear estimates 

In this section, attention is now given to L^ — ?■ L°° estimate of solutions to the linear 
equation of (11. 2p . 



Theorem A.l. Consider the linearized equation of (11. 2p 
(A.l) Ut-Au = 

with u{0) = uo= 1^°]. IfuoeX^n{L^x L^), then u{t) G X^ fl (L°° x L°°) and satisfies 



(A.2) 



\u{t)\^^oo < (1+i) Kll^olll + |%|l,l), 



for t > 0. 



To prove Theorem lA.il we need the following two lemmas. The first one is the well-known 
Van der Corput lemma which we state without proof [291 ED] • 

Lemma A.l. Let h be either convex or concave on [a, b] with — oo < a < b < oo. Then 



^h(Od^ 



< 4{min \h"\}-^/^ for h" ^^ in [a, b]. 

[a,b] 



Applying the above Van der Corput lemma to the linearized equation f lA.ip . we obtain 
Lemma A.2. Let nA > 1 and < e < 1. Then 



(A.3) 

where h{^, a) 



sup 

a6R 



Jth(i,a) 



d^ 



I5I<" 



< (e + t 2 niax{e 2,722} 



e 



+ <. 



Proof. We use the stationary phase method [261 [29] and Lemma lA.ll Since 

dh 1 - e^ 



d^h 2e(e' - 3) 



and 



de 



:i+e 



2^3 
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for e small and n large enough, it follows from Lemma 2.2 that 



,ith{i,a) 



€<|5|<V3-£ 



di 



< 4t 2 min 






e,a 



d^h 



de 



(73 



and 



/ 



Jth{S,,a) 



VS+e<\^\<n 



On the other hand, 



and 



d^ 



< 4t 2 mill 



d'^h 



3 1 



ra, a, 



d'^h 



e, a 



(73- 



< ce 2t 2 



5^2V- ' - 



e, a 



< ct 2 max{n2,e 2|. 



l?l<6 



< ce 



,jth(?,«) 



'V3-e<|5|<V3+e 

Therefore, the above estimates imply flA.31) . 



di 



< ce. 



Proof of Theorem \A.l\ One can write the solution u = S{t)uo of (lA.ip as 
(A.4) uit) 

where Uq is the Fourier transform of Uq. It follows that 






|M(i)|oo,oo < CS 



(r/0 ± ^o)e^*''^^'^"«/*)rfe 



where the sum E is over all two sign combinations. Using Lemma IA.21 we obtain 



\u{t)\ 



00,00 r^ 



(l^ol + \vo\)di 



\i\>n 



(r/o±^o)e^*'^«'^^)rfe 



ICK" 



^ ll%lli 



[i+evdi 



'\i\>n 

+ / (l^o(a;)| + \vo{x)\)dx 



jthi(,±^) 



d^ 



|C|<n 



^ -ill-, II / _i r -i i^ T \ 1^ I 

^n 2 lluolli + ( e + t 2 max I e 2,71,2} 1 |no|i,i. 



D 



Therefore the estimate (]A.2p can be obtained by choosing e = t 3 and n = t*. Since S'(t) 
is a unitary group in X* for any s > 0, we know u{t) = S{t)uo G X^. D 
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